Soliton mode- locking in optical microresonators 
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The discovery of passive mode-locking via saturable absorbers [Tl [2] has led to optical 
femto-second pulses [3, 4] with applications ranging from eye-surgery [5] to the analysis 
of chemical reactions on ultra-short timescales[6j. In the frequency domain a train of 
such optical pulses corresponds to a frequency comb (equidistant optical laser lines 
spaced by the pulse repetition rate) [7— 9], which find use in precision spectroscopy and 
optical frequency metrology [ ToHl2] . Not relying on mode- locking, frequency combs 
can also be generated in continuously driven, high-quality-factor Kerr-nonlinear op- 
tical microresonators via cascaded four- wave mixing [13] . Over the past years these 
Kerr-combs have been demonstrated in a variety of microresonator geometries |14H20| . 
Applying a pulse-shaping mode-locking mechanism could enable compact and robust 
femto-second pulse generators. However, conventional and emerging (e.g. Graphene 
based) saturable absorbers are challenging to apply to microresonators, as they af- 
fect the high-quality-factor essential to nonlinear frequency conversion. Here, we re- 
port on passive mode-locking in microresonators without saturable absorber. This 
mode-locking is achieved via soliton formation [Tl [2] supported by the balance between 
anomalous resonator dispersion and Kerr-nonlinearity induced self-phase modulation. 
The transitions to and between different soliton states manifest themselves itself as dis- 
crete steps in the resonator transmission. Numerical modeling of the nonlinear coupled 
mode equations is used to physically understand the mode-locking and in combination 
with an analytical description allows the identification of mode-locked regimes. Exper- 
imentally, we observe the generation of pulses with 200 fs duration. Equally important 
low noise and low line-to-line power variation comb spectra are achieved via single soli- 
ton states. These spectra bring microresonator combs decisively closer to frequency 
domain applications such as channel generators in advanced telecommunication or in 
fundamental studies such as astrophysical spectrometer calibration[13j. In the time 
domain the presented results open the route towards compact femto-second sources, 
where the broadband parametric gain in principle allows for sub-cycle pulses. More- 
over, femto-second pulses in conjunction with external broadening provide a viable 
route to microresonator RF-to-optical links[9l 121] . 



Frequency comb generation in optical microresonators [22 is achieved via cascaded- four- wave mixing (FWM, cf. 
Fig lc) mediated by the Kerr-nonlinearity [14]. Here, a continuous- wave (CW) pump laser is converted into equally 
spaced optical modes, where the mode spacing corresponds to the free-spectral range (FSR), or equivalently, the inverse 
resonator roundtrip time 1/Tr. Driven by a wide range of application in telecommunications [23 , spectroscopy [24]. 
astronomy[25, 26 , and compact space optical clocks[27 this new class of Kerr-combs has rapidly evolved in various 
materials and geometries [28 . Significant recent advances include the demonstration of octave spanning spectra [29] [30], 
full phase stabilization [19] [31], low phase noise microwave generation [15] [19] [31] [32], fully planar CMOS compatible 
design [T6 ] fT7 ] [33] . extension of spectral coverage from visible [34 to mid-infrared wavelength [35 , as well as arbitrary 
optical waveform generation [36]. The parametric four- wave mixing process results in coupled phase relations between 
all the optical modes, which leads to the generation of a periodic time domain output (cf. Fig Id). Importantly 
however the output does not correspond intrinsically to pulses. External line-by-line phase and intensity adjustment 
have been used after comb generation for pulse-shaping [32] [36], but is restricted to a small number of comb modes. 
Mode-locking, i.e. full phase synchronization inside the microresonators would naturally yield a train of optical 
pulses with a pulse repetition rate f rep = 1/Tr (cf. Fig Id) and enable ultra-short pulse generation for widespread 
application. Moreover, mode- locking would provide a way to achieve low- noise frequency comb spectra, with little 
line-to-line variation and without spectral gaps - a goal difficult to reach in systems with mode-spacings below 100 
GHz, where multiple, inconsistent subcombs may form [37] . 

Introducing a saturable absorber or an equivalent optical element into the laser cavity [TJ [2], easily possible in 
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conventional laser systems is hindered by the requirement of a high optical quality factor Q. Still, mode- locking 
mechanisms intrinsic to microresonators could be present. For instance self-focusing [TJ [2], affecting both the intermode 
coupling inside the resonator as well as the input output coupling to the waveguide. Self- focusing could also induce 
an intensity dependent surface loss. Moreover, a small saturable component of material absorption could facilitate 
mode-locking. Interestingly, pump laser detuning dependent transition from high to low phase noise comb states 
have been observed [20j |30l [37j [38] and interpreted as a signature of potential phase locking [20j |30l [38], but an 
understanding is still lacking. Analysis of the damped, driven non-linear Schrodinger equation (NLS) suggests that 
mode-locked regimes could exist in microresonators [39]. It has also been proposed that Kerr-cavity solitons may form 
in microresonators [40 J. This soliton formation could lead to mo de-lo eking [TJ [2]. However, to date no definitive proof 
nor mechanism of mode-locking has been identified. 

To reveal potentially mode-locked states as a function of laser detuning we scan the pump laser uj p (in the direction 
of decreasing optical frequency) over a high-Q resonance (Q=400 million, free-spectral range FSR=35.2 GHz) of a 
crystalline MgF 2 resonator [35j I37 J |4T | l42](cf. Fig. la). The radio frequency (RF) signal generated by the beating 
between neighboring comb lines (down-mixed to 20 MHz) serves as an indicator of mode-locked states. Practically, we 
sample and subsequently Fourier-transform the RF signal, while the laser scan is performed. In mode-locked states 
the resulting spectral beat signal is expected to be a single low-noise peak. 

Fig. 2b shows the evolution of the optical spectrum during the laser scan. Reducing the laser-cavity detuning leads to 
a build-up of intra-cavity power and once a critical power threshold is reached widely spaced primary comb modes are 
generated, followed by secondary modes filling in the spectral gaps as frequently observed (cf. Fig. lc)[36, 37, 42 . The 
increase of intra-cavity power also results in a shift of resonance frequency towards lower frequencies via a combined 
effect of thermal expansion, thermal refractive index change and Kerr-effect. This leads to a non-Lorentzian, triangular 
resonance shape when the pump laser is tuned with decreasing optical frequency over the resonance (thermal triangle, 
cf. Fig 2a, inset) [43, 44 . Importantly, the thermal effect enables self- locking of the cavity to the laser [44 . 

The Fourier-transformed, sampled RF signal is contained in Fig. 2c. We observe a transition from a broad, noisy 
RF signal to a single, low- noise RF beat note. Strikingly, this transition coincides with the beginning of a series of 
discrete steps in the transmission, which deviate markedly from the expected thermal triangle. Such series of discrete 
steps have been observed in non-linear microwave resonator, and connected to soliton for mat ion [45]. These 
observations are highly indicative for mode-locking. 

To investigate the intriguing observations of discrete steps in the transmission, and the possible connection to mode- 
locking we carried out numerical simulations based on the coupled mode-equations approach recently introduced to the 
field of microresonator combs [46]. The full control over the included physical effects allows identifying the underlying 
physical mechanisms. Here the simulated system corresponds to a typical MgF 2 microresonator with resonance 
frequencies uj^ = ujo + /i • D\ -f /i 2 • D2/2 -f/i 3 • ^3/6, where D\ = 2tt • FSR, D2 and D3 correspond to the FSR, and the 
second and third order dispersion around the resonance frequency ujq of the pumped mode [37 having the mode index 
H = 0. The remaining resonator simulation parameters are refractive and non-linear indices, as well as, the effective 
mode- volume (cf. methods). We neglect effects of non- unity mode-overlap, interactions with other mode families, 
any particularities of the resonator geometry and thermal effects. The resulting set of coupled mode equations (cf. 
methods) is propagated in time using an adaptive step-size Runge-Kutta integrator. Results of a numerical simulation 
including 101 optical modes are shown in Fig. 3a,b,c. The blue curve in panel a shows the simulated intracavity 
power as function of ( Q = 2(uo — uo p ) / n , which is a normalized measure of the laser tuning, where k = uoq/Q denotes 
the cavity decay rate. The blue curve is mirrored horizontally for convenience when comparing to the equivalent 
experimental transmission trace in Fig. 2a (an increased transmission corresponds to a drop in intracavity power). 
Remarkably, the step features are very well reproduced, implying that the simulation fully includes the relevant 
physical mechanism. In agreement with the experiment the number and height of steps fluctuate in repeated scans. 
This last part of comb evolution is sensitive to fluctuations in the pump power and, numerically, short pump power 
drops can be used to induce a transition to a new comb state. Numerically tracing out all possible comb evolutions 
yields the orange curves in Fig. 3a. The evolution of the optical spectrum during the scan in Fig. 3a is shown in Fig. 
3b and follows the typical path [37]. With each steps the optical spectrum becomes less modulated until it eventually 
reaches a perfectly smooth envelope state (frame XI). 

To reveal mode-locking we investigate the time dependent waveform in Fig. 3c by summing the simulated individual 
optical modes. Indeed, the first step (frame V) corresponds to a transition to a state where multiple pulses inside the 
cavity exist. Further steps can be associated with a stepwise reduction of the number of pulses propagating in the 
resonator. The separation between pulses is random. The presence of anomalous dispersion and Kerr-nonlinearity in 
the cavity (cf. Fig. lb) suggests that the pulses may actually be bright, temporal cavity solitons. These have been 
observed for the first time in a continuously driven fiber-loop cavity and proposed for the case of microresonators 
frequency combs[40]. To confirm the soliton nature of these pulses we perform a simulation of 501 modes (cf. Fig. 
3d,e,f) and analyze a state of five pulses. We compare the numerical simulation with an approximate analytical 
solution for the damped, driven non-linear Schrodinger (NLS) equation given by a soliton on a CW background [47] 
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(Note that this NLS equation is equivalent to the coupled mode- equations when neglecting third order dispersion 
[39 J. Based on this, we derive the following expression for multiple solitons 

* w -I + + 'V 2< °-5?) f>» (7I<* - *>) ■ w 

where -0 denotes the complex field amplitude, the angular coordinate inside the resonator, <pj the angular coordinate 
of the jth soliton, u p the pump power, N is the number of solitons, and / corresponds to the pump power (cf. methods). 
We note that this solution, except for the soliton positions ^, does not include any free parameters. Indeed, the close 
to perfect match between analytical solution and numerical result shows, that the pulses forming in the microresonator 
are temporal ca vity soli tons. Eq. [I] allows to estimate the detuning and dispersion dependent temporal width of the 
soliton T s = y // ^2/2Co/^i- Note that due to the non-linear Kerr-shift of the resonance frequency, the effective 
detuning between pump laser and resonance is smaller than (o- 

Having revealed the soliton nature of the pulses in numerical simulation, we can interpret the blue curve in Fig. 
3a based on general limits applying to solitons as solutions of the driven, damped NLS [48 . Adopting these criteria 
for the present case we identify three main regions colored red, yellow and green in Fig. 3a. Solitons with a constant 
temporal envelope can only exist in the green area. While the yellow area still allows for solitons with time varying 
envelope (breather solitons) [49] , solitons can not exist in the red area. The curves defining the yellow area are given 
by |a| 2 ± = | Co ± | vCo 3, where \a\ 2 refers to the average energy in the resonator. An additional restriction limits 
detuning to Co = K 2 /8f 2 , beyond which no solitons can exist [48] . Note that in the red area on the left the system 
may undergo chaotic Hopf-bifurcations[48 . The correspondence of these analytical limits with the numerical results 
is remarkable. 

By averaging the analytical soliton solution (eq. [I]), for different number of solitons, over one cavity roundtrip time 
we can derive the total internal field inside the resonator. The result is shown as dark gray dashed lines, which are in 
excellent agreement with the numerically observed steps (to account for the limitation due to the low mode number 
an additional correction factor of order unity is applied). The step height is directly proportional to the number of 
solitons removed in a comb transition. 

The simulations show that microresonators can intrinsically be mode- locked via soliton- formation [1 , 2 , which relies 
on the anomalous cavity dispersion and the Kerr-non-linearity causing frequency and intensity dependent phase shifts. 
Other physical mechanism, particularly an additional saturable absorber, are not included in the simulation. For 
completeness we note that the soliton transition may also be observed in simulations of an injection locked operation 
scheme, where the pump power is gradually increased and the laser remains at zero detuning with respect to the 
non-linear ly shifted resonance. 

Having reached a theoretical understanding of the soliton states, we next experimentally investigate - in addition 
to their RF beatnote and optical spectrum - their temporal characteristics by performing a frequency-resolved optical 
gating experiment (FROG, Fig. 4a and 4e for the setup) [50]. This corresponds to a second-harmonic generation 
(SHG) autocorrelation experiment, where the frequency-doubled light is spectrally resolved (cf. methods). In full 
consistency with the numerical simulations, we observe single and multiple pulse states as shown in Fig. 3c. The 
one pulse state is characterized by a smooth spectral envelope, without spectral gaps. The power spectral envelope 
exhibits a squared hyperbolic secant sech 2 -shape (with FWHM = 1.6 THz) as expected from the Fourier transform 
of a sech-shaped soliton pulse. Based on the time-bandwidth-product of 0.315 the expected pulse duration is 197 fs. 
The low phase noise RF beatnote is resolution bandwidth limited to 1 kHz and its signal-to- noise ratio exceeds 60 dB. 
The FROG trace shows pulses well separated by the cavity roundtrip time of Tr =28.4 ps, corresponding to the FSR 
of 35.2 GHz. The multi-pulse states (here shown for the case of two and five pulses), show a more structured optical 
spectrum. This structure reflects the number and distribution of pulses in the cavity. The RF beat note generated in 
the multi-pulse states is of similar quality as in the single pulse state. Importantly, the FROG measurement allows 
for a full reconstruction (neglecting a time direction ambiguity) of intensity and phase of the pulses (cf. Fig. 4b). 
The reconstructed intensity is consistent with the expected sech 2 -shape for solitons and has a FWHM of 200 fs in 
agreement with the expectation. The FROG traces show that it is the full optical spectrum that contributes to the 
pulses. This is particularly true in the case of multiple pulses. 

Combining experimental, numerical and analytical results, we have demonstrated mode-locking in a MgF2 microres- 
onator via soliton formation, which is only based on the generic properties of Kerr nonlinear microresonators with 
anomalous dispersion. Consequently, our observations are generic and apply equally to other Kerr comb platforms. 
From a frequency domain perspective, soliton mode-locking enables microresonator based frequency combs with low 
noise and smooth spectral envelope. Such low line-to-line variations are key to applications such as channel generation 
for advanced modulation formats in telecommunication. From a time domain perspective it enables ultra-short soliton 
pulse generation in a microresonator. In combination with chip-scale [16] [17] integration this opens the route towards 
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compact, low cost ultra-short pulsed lasers sources . Due to the broadband nature of the parametric gain, which can 
significantly exceed the spectral width of conventional laser gain media, the direct synthesis of sub-cycle pulses with 
this method is conceivable. 



I. METHODS 

Kerr- Comb generation: The pump laser used is a narrow linewidth (100 kHz shortterm) 1553 nm fiber laser 
amplified by an erbium doped fiber amplifer (EDFA). Evanescent coupling of the resonator is achieved via a tapered 
optical fiber [51] yielding a coupled linewidth of nj2i\ = 450 kHz. 

Numerical simulation: Besides the resonance frequencies uj^ defined by Di/2tt = 35.2 GHz, D2/2tt = 10 kHz, 
Ds/2tt = —130 Hz as measured [52 , the resonator is characterized by its refractive and non-linear Kerr indices n = 1.37 
and n 2 = 0.9 x 10" 20 m 2 W _1 , the effective mode volume^ V eff = 5.6 x 10" 13 m 3 and quality factor Q = 250 x 10 6 
. Starting from the coupled mode equations for the non-equidistant cavity modes at the resonance freqeuncies uj^ 
a set of dimensionless coupled mode equations for the equidistant comb modes a M at a priori unknown frequencies can 
be derived (cf. SI in [37]). The mode amplitudes are related via the transformation a M = A fI ^2g/ K,e' l< ' UJ ^~ UJp ~ flDl ^ t 
where |A M | 2 is the number of photons. The coupled mode equations read: 



fj/ — — K /x" >/V 



The Kerr non-linearity is contained in the nonlinear coupling coefficient g = ^2y™ 2 f •> f = y is the nor- 

malized driving field, and r = nt/2 the cavity intrinsic time scale. 2K + 1 is the number of simulated modes 
and C/i = 2(a; M — u p — fiDi)/ n. In a post-processing step the optical frequencies, field amplitudes, phases and RF 
signal can be inferred by taking amplitude and phase modulation of the simulated field amplitudes a M into account. 
Random vacuum fluctuations are used to seed initial FWM. The pump power is set to P in — lOOmW at a pump 
frequency of u p = 193 THz. The computation time scales cubically to highest order and a simulation of 501 modes 
takes 5 days on a standard PC. 

FROG Experiment: Prior to the FROG experiment [50 the optical spectra are sent through a fiber-Bragg grating 
for pump supression (30 dB) and subsequently amplified to 50 mW. Dispersion compensating fiber (DCF3,DCF38) 
is used for approximate dispersion compensation. In the FROG setup (cf. Fig 4c) the generated optical pulses, are 
interferometrically split and recombined with an adjustable delay in a non-linear crystal resulting in self-gated second 
harmonic generation (SHG), that is doubling of the optical frequency whenever the optical pulses in the two arms of 
the interferometer overlap temporally in the non-linear BBO crystal. The generated SHG light is spectrally resolved 
and recorded as a function of delay, yielding a so called FROG trace. Each scan consist of nearly 1000 spectra with 
individual exposure times of 800 ms. The higher resolution N-by-N (N=63) FROG trace of the single pulse state in 
Fig. 4b) is analyzed using a principal component generalized projection FROG algorithm, after noise removal via 

Fourier-filtering. The FROG reconstruction error is defined as e = y^ X^fj (M™ eas — M^ eco ) 2 , where M™ eas and 
M^ co denote the elements of the N x N matrices representing the measured and reconstructed FROG traces. 
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FIG. 1: Fig.l Mode-locking and frequency comb generation in microresonators. a. MgF 2 crystal carrying two 
whispering-gallery- mode microresonator of different size (here: the smaller one with an FSR of 35.2 GHz is used). Here an 
optical whispering-gallery-mode propagates along the circumference L of the resonator within the roundtrip time Tr. The 
magnified panels show a zoom in and the simulated optical mode profile, b. Anomalous dispersion and Kerr-effect leading 
to a frequency uj dependent and intensity / dependent effective refractive index n e ff (no: refractive index; 712: nonlinear 
refractive index), c. Frequency comb generation in Kerr-nonlinear microresonator via four- wave- mixing, The newly generated 
optical lines are supported by the cavity's resonances (gray). Widely spaced primary lines emerge first, followed by secondary 
comb lines filling the gaps. d. Optical whispering gallery comb modes (blue) and resulting intensity (red). If the relative 
phases (indicated by the black arrows) are constant but random the resulting intensity is periodic with Tr; if the phases are 
synchronized via mode-locking the resulting intensity is pulsed with a pulse repetition rate of 1/Tr. 
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FIG. 2: Fig. 2 Transmission and beat note a. Transmission observed when scanning a laser over a high-Q Kerr-non- linear 
resonance in a MgF 2 resonator (coupled pump power 5 mW). The transmission signal follows the expected thermal triangle 
(cf. inset) with deviations in the form of discrete steps (green shading), b. Evolution of the optical spectrum for three different 
positions in the scan (I, II, III), c. Downmixed radio frequency (RF) beat signal, d. Experimental setup composed of generic 
coupling setup (including pump laser and resonator) followed by an optical spectrum analyzer (OS A), an oscilloscope to record 
the transmission and to sample the down-mixed beatnote (via the third harmonic of a local oscillator LO at 11.7 GHz), and an 
electrical spectrum analyzer (ESA) to monitor the beatnote. Before beatnote detection the pump is filtered out by a narrow 
fiber-Bragg grating (FBG) in transmission (FPC: Fiber polarization controller; CW continuous wave; EDFA: erbium-doped 
fiber amplifier; PD: photodiode). 
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FIG. 3: Fig. 3 Numerical Simulations a. Transmission (blue) in a simulated laser scan over a resonance in a MgF2 resonator 
leading to step-like features (101 simulated modes). The orange lines trace out all possible transmission evolutions of the system. 
The dashed lines show an analytical description of the steps. The green area corresponds to the area, where solitons can exist, 
the yellow area allows for solitons with time variable envelope; no solitons can exist in the red area. b/c. Optical spectrum and 
intracavity intensity for different positions I-XI in the laser scan. d. Optical spectrum obtained when simulating 501 modes 
and stopping the laser scan in the soliton-regime. e. Intracavity intensity for the comb state in (d) showing 5 solitons. f. 
Zoom into one of the soliton states showing the numerical results for the field real (red) and imaginary part (dark blue). The 
respective analytical soliton solutions are shown in dashed light blue and orange. 
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FIG. 4: Fig. 4 Experimental demonstration of soliton states a. Optical spectra (left column) of three select comb states 
(one, two and five pulses). The insets show the RF beat note, which is resolution bandwidth limited to 1 kHz width in all cases. 
The FROG (frequency resolved optical gating) traces (right column) reveal single and multiple pulses per cavity roundtrip time. 
The red line in the optical spectrum of the one pulse state shows the spectral hyperbolic-secant envelope expected for soliton 
pulses with a full- width- half-maximum (FWHM) of 1.6 THz. b. Higher resolution experimental FROG trace of a one soliton 
pulse (left). The reconstruction converges to a FROG error of 1.7% in good agreement (middle) with the experimental trace. 
The reconstruction yields an estimated pulse duration of 200 fs (FWHM). c. Setup for the FROG experiment: After pulse 
generation in the microresonator (cf. Coupling Setup, Fig. 2) the pump light is filtered out using a fiber-Bragg grating (FBG) 
in transmission. The optical signal is dispersion compensated via dispersion compensating fiber (DCF3, DCF38) and amplified 
in an erbium doped fiber amplifier (EDFA) to approx. 50 mW. The amplified light is free-space coupled to an Michelson-type 
interferometer with offset beams. Focused by a lens the two beams overlap in a non-linear BBO crystal for frequency doubling. 
Only the frequency doubled beam in the central forward direction is allowed to reach a grating based CCD-spectrometer for 
detection. The delay of the interferometer is varied to record the frequency doubled optical signal spectrally resolved as a 
function of delay (HWP: half- wave-plate, QWP: quarter- wave-plate) . 



